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Electron emission from hafnium carbide (HfC) field emission tips induced by a sub- 10 fs, 150 MHz 
repetition rate Ti: sapphire laser is studied. Two-photon emission is observed at low power with a 
moderate electric bias field applied to the tips. As the bias field and/or laser power is increased, the 
average current becomes dominated by thermally-enhanced field emission due to laser heating: both 
the low thermal conductivity of HfC and the laser's high repetition rate can lead to a temperature 
rise of several hundred Kelvin at the tip apex. The contribution of current from a thermal transient 
at times shorter than the electron-phonon coupling time is considered in the context of the two- 
temperature model. Under the conditions of this experiment, the integrated current from the 
thermal transient is shown to be negligible in comparison with the two-photon emission. A finite 
element model of the laser heating and thermal conduction supports these conclusions and is also 
used to compare the nature of thermal effects in HfC, tungsten, and gold tips. 

PACS numbers: 79.60.Jv, 78.47.J-, 78.67.-n, 81.70.Pg 



I. INTRODUCTION 

Field emission tips are extremely high brightness 
emitters^ with large transverse coherence length^ and 
small energy spreacP. They are therefore excellent 
sources for interferometry experiments such as elec- 
tron holograph};^ and coherent electron diffraction.^ 
Femtosecond field emission sources are being devel- 
oped by triggering electron emission with ultrafast laser 
pulses .^^H^ These sources have extremely small emis- 
sion areas in co mparison with existing ultrafast elec- 
tron sourceJ^^' based on flat photocathodes and should 
therefore lead to better resolution for diffraction and 
imaging. Furthermore, the small source size leads to the 
high transverse coherence required for electron interfer- 
ometry and holography.!^ Ultrafast laser-triggered tips 
promise to extend the full range of electron microscopy 
imaging techniques to the femtosecond time-scale, creat- 
ing tools for the study of chemical and physical transfor- 
mations of materials with atomic-scale resolution. 

A number of areas have already benefited from laser- 
triggered emission from sharp tips. Optical field enhance- 
ment at the tip apex has allowed exploration of strong 
field physics in solid state systems without the need for 
laser amplifier s.^1^ Due to the extreme sensitivity of 
current to the light intensity at the emission site, field 
emission tips have been used to study the temporal dj^ 
namics of the optical near- field induced by the tip itself^ ^ 
and map optical fields in nearby nano-structures^ with 
nanometric spatial resolution. Ultrafast laser-triggered 
emission has also been used as a source of electron 
pulses for experiments that verified the dispersion-less 
nature of the Aharonov-Bohm effect and is being devel- 
oped as an electron source for ultrafast scanning electron 
microscopy. These sources could also enable laser-based 
electron accelerator^ and a table-top ultrafast x-ray mi- 
crofocus source. 



Hafnium carbide (HfC) is an attractive material to use 
for ultrafast laser-induced electron emission because of 
its low work function (~ 3.4 eV), which makes it easier 
to liberate electrons, and its good mechanical stability at 
high temperatures.^^ The low work function relaxes laser 
intensity requirements for electron emission, which would 
allow devices that use low power fiber lasers to excite the 
electrons or devices where electron optics elements near 
the tip preclude high numerical aperture focusing of the 
laser. Here, we confirm this by showing that to generate 
one electron per pulse (an important regime for avoiding 
space charge while maintaining the maximum possible 
flux) HfC requires one quarter of the pulse energy as re- 
sults reported in tungsten^^ and gold.^^ HfC has been 
demonstrated to have higher emission current densities 
than tungsten in stable DC operation,"^ and it may be 
annealed more aggressively than tungsten to create clean 
surfaces. The low thermal conductivity of HfC is, how- 
ever, a disadvantage if higher repetition frequencies or 
greater numbers of electrons per pulse are desired; it re- 
sults in higher average temperatures and longer-lasting 
elevated temperature transients for pulsed laser heating 
in comparison with materials like gold or tungsten. 

Several ultrafast laser-induced electron emission pro- 
cesses have been observed in field emission tips, and these 
in general transfer the precise timing of the laser pulses to 
the electrons. In the perturbative regime, photo-assisted 
field emission^ and multi-photon emission^ are prompt 
with res pect to the laser intensity envelope. Work in 
tungsteipl^ and gold^^ has explored higher pulse flu- 
ences, where for sufficiently short optical pulses the 
emissio n pro file is predicted to have sub-optical cycle 
features.L^^ 

Here, we show that laser-induced electron emission 
from HfC is due to two-photon emission over a broad 
range of laser power and bias field, yielding up to several 
electrons per pulse at 150 MHz. However, for certain pa- 
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rameters we observe thermally-enhanced field emission 
(TFE), which we describe as a transient current pulse 
(transient TFE) superimposed on a steady-state back- 
ground (DC TFE). DC TFE arises from a significantly el- 
evated average apex temperature due to the low thermal 
conductivity of HfC and the laser's high repetition rate. 
A transient TFE pulse could arise from rapid heating of 
the electrons by each laser pulse, followed by electron- 
phonon coupling on a ps timescale and cooling by con- 
duction. The sub-ns timescale of conduction following 
laser-induced heating is known to play a role in ultra- 
fast laser- assisted field evaporation,^^ which is used for 
atom-probe tomography. However, until now there have 
been no investigations of the effects of a laser-induced el- 
evated electron temperature on the timing properties of 
electron emission from fs-laser-triggered tips. By com- 
paring the scaling of DC and transient TFE with that 
of multiphoton emission, we show that over the parame- 
ter range of our experiment, the transient effect is not a 
dominant source of current. 

The paper is organized as follows. The experimen- 
tal setup is described in Section [TTj and Section [TTT| is an 
overview of emission processes. Section IV A presents ev- 
idence for two-photon emission and DC thermal-assisted 
field emission as the dominant mechanisms observed in 
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HfC. Section IV B describes modeling of laser heating 
of the tip. In section IV C[ these results are combined 
to place an upper bound on the fraction of pulsed cur- 
rent that is due to transient TFE and to compare how 
multiphoton emission and transient thermal-assisted field 
emission scale with pulse energy in a low repetition rate 
system. 



II. EXPERIMENTAL SETUP 

A diagram of the experimental setup is shown in Fig. 
[l] The output of a sub-10 fs, 150 MHz repetition fre- 
quency Tiisapphire oscillator is focused with a 2.7 pm 
waist (1/e^ intensity radius) onto a 310-oriented hafnium 
carbide field emission tip (Applied Physics Technologies, 
McMinnville, OR). The dispersion of the vacuum win- 
dow and other glass in the beam-line is compensated with 
double-chirped mirrors, and the ~ 10 fs pulse duration 
at the tip is verified by an interferometric autocorrela- 
tion with electron emission from the tip as the non-linear 
detector. The tip has a 120 nm radius of curvature at 
the apex. We take the work function of the tip's emitting 
surface to be 3.4 eV.^^ The tip is maintained in ultra-high 
vacuum (< 5 x 10~^^ torr). Surface contamination modi- 
fies the work function of the tip; we therefore flash anneal 
to 2100 K for 1 s prior to each set of measurements. A 
DC bias is applied between the tip and an extraction 
aperture, and the spatial distribution of the electrons is 
imaged using a dual chevron configuration microchan- 
nel plate (MCP) with a phosphor screen. The absolute 
gain of the MCP was determined by counting detection 
events in high gain operation and then using a relative 




FIG. 1. Experimental setup. The ultrafast laser pulses are 
focused onto a HfC tip, which is maintained in ultra-high vac- 
uum (< 5x 10"^*^ torr). A bias voltage, Vtip is applied between 
the tip and extraction aperture, located 1.5 mm away. On the 
other side of the aperture, the electrons propagate in a field 
free region. Electrons are detected with an MCP; a bias tee 
splits the RF and near DC components of the current at the 
front of the MCP. A chopper in the laser beam allows lock-in 
detection of luser- 



calibration for lower voltages across the plates. 

A bias tee splits the RF and near-DC components of 
the current at the front (incident surface) of the MCP. 
The DC component is used to measure the current in- 
duced by the laser (/laser) with lock-in detection at the 
frequency of a chopper in the optical beam path, /laser is 
equivalent to the difference between the average current 
incident upon the MCP with and without light on the tip. 
We use /dc to describe current from the tip when there 
is no laser illumination. The RF signal contains harmon- 
ics of /rep, the laser repetition frequency, extending to 
frequencies > 2.5 GHz, limited by the MCP's response. 
Since the electron detection events have a known pulse 
shape, the amplitude of the harmonic at 2/rep can be 
converted into an equivalent average laser-induced cur- 
rent which we define as /2f. For instance, when all the 
signal is pulsed (no appreciable DC TFE), /2f and /laser 
have the same value. 



III. EMISSION PROCESSES 

A diagram of DC and laser-induced electron emission 
processes is shown in Fig. |2j We summarize each of these 
below. 



A. Field Emission 

The field at the apex of a metal tip is given hy F = 
V/(/cr), where r is the radius of curvature at the apex of 
the tip and A: is a geometrical factor (typically ~4-5). For 
tips with r of 10-100 nm, field strengths of GV/m may be 
obtained using modest voltages. Under such large bias 
fields, electrons tunnel directly from the fermi level into 
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FIG. 2. Electron emission processes, (a) Field emission (in- 
cluding thermally-enhanced), (b) Photo-assisted field emis- 
sion, (c) Multi-photon (over the barrier) emission, (d) Above- 
threshold photoemission and optical field emission. 



the vacuum, as shown in Fig. |2|a). Without illumina- 
tion, the magnitude of the total current density emitted 
from a field emission tip as a function of field, F, and 
electron temperature, Tei, is given by 



J(F,Tei) = e j D{E,,F)N{E,,T,i)dE,, 



(1) 



where e is the elementary charge. Ex is the energy of 
the electron in the surface normal direction, D{Ex^F) 
is the transmission coefficient of the tunnel barrier, and 
N{Ex^Tq\) gives the distribution of electrons with nor- 
mal energy Ex at a given temperature. The problem is 
analyzed with x as the surface normal direction. The po- 
tential inside the metal is approximated as flat, and for 
X > 0, the potential energy barrier seen by the electron 
is simply 



U{x) = (j) ^ ji — eF X ■ 
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(2) 



referenced to zero inside the metal, with (j) the work func- 
tion of the metal, and fi the fermi energy. The final term 
is due to the image potential. Using the WKB approxi- 
mation, 



D{E,,F) ^ exp (-2 VW) 



Ex dx 



where x+ and x_ are the classical turning points. We 
take D{Ex^F) ^ 1 when Ex exceeds the barrier height. 
For a free electron metal. 



N{Ex,T,,)dEx 



\n{l^e^^-^-^^^^^^')dEx, 
(4) 

with h the Planck constant and the Boltzmann con- 
stant. This simple model for tunneling is extremely suc- 
cessful in analyzing the current vs. voltage {I-V) charac- 
teristics of field emitters at room temperature.^ At zero 
temperature, the integral can be carried out analytically 
to obtain the Fowler-Nordheim equation: 



,3^2 



S7Th(l)t^{w) 



exp 



3heF 



-v{w) 



(5) 
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Here, w = e^/^ Y^F/47reo/^, and t{w) and v{w) are rel- 
atively slowly varying functions of order 1, related to 
the elliptic integrals that occur in Eq. ([l]).^ Using this 
in conjunction with Eq. ([5| multiplied by the emitting 
area, one obtains a relationship between the measured 
quantities / and V. Fig. |3|b) shows a Fowler-Nordheim 
equation fit to room-temperature I-V data from our tip 
(no laser illumination) along with the corresponding fam- 
ily of I-V curves (/tfe(F, F)) at elevated temperatures, 
which have been obtained by numerically integrating Eq. 
([T]). Increasing the temperature populates states above 
the Fermi level which enjoy an exponential increase in 
the tunneling transmission, resulting in enhanced emis- 
sion that can be very sensitive to temperature as well as 
field. We refer to emission at elevated temperatures (so 
that the Fowler-Nordheim equation is no longer a good 
description) as thermally-enhanced field emission (TEE). 
Note that, as we use it, this term includes thermionic 
emission.^ For the parameters of our experiment, the 
temperature nonlinearity can scale locally as steeply as 

2^12 



B. Ultrafast laser-induced emission 

With a mode- locked laser, it is possible to achieve both 
high peak intensities and few cycle pulse durations. This 
light, incident upon a biased field emission tip, leads to 
enhanced electron emission that can be described by the 
processes m Fig. ^h-d). Proper management of these 
processes should allow the temporal profile of the pulse 
train, or even of each optical cycle, to be transferred to 
the temporal profile of electron emission. 



1. Multiphoton emission 

For visible hght, the energy of a single photon is typ- 
ically less than the work function of the metal. Absorp- 
tion of multiple photons can liberate electrons over the 
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barrier (Fig. |2jc)). A''-photon emission is described by 



Jn =ANr [{NTwj - 0)^ + y (fcsTei)^ 



(-1) 



n+1 



-n{Nnoj-(i))/kBTe\] 



(6) 



n=l 



with /, laser intensity, cj, angular frequency of the laser, 
and An a constant. Multiphoton emission has been ob- 
served in gold^ and tungsten^^ tips, with N = A and 
A/" = 3, respectively (both experiments used a Tiisapphire 
laser). When a bias field is applied to the tip, the 
height of the potential barrier is reduced due to the im- 
age charge when an electron leaves the metal (Schottky 
effect).^ To account for this effect, we make the substitu- 
tion (j) ^ (j)- Y^eF/(47reo) in Eq. [g] 



where m, e are the mass and charge of the electron, re- 
spectively. The photo-field and multiphoton emission 
regimes occur for 7 ^ 1. The quasi-static or 7 <C 1 
regime is interesting from the perspective of ultrafast 
electron sources because, since tunneling can only hap- 
pen during one half of the optical cycle and is extremely 
nonlinear in field, electron pulses that are shorter than 
the optical cycle of the laser may be generated. Experi- 
ments with tungsten tips have reached 7 ~ 1 which 
is possible with 10-100 mW of laser power at ~ 100 MHz 
repetition rates due to optical field enhancement by the 
tips. Even in this "intermediate" regime, numerical sim- 
ulation of the Schrodinger equation predicts sub-optical 
cycle features in the electron emission.l^ For atoms, this 
has been predicted analytically by an extension of the 
Keldysh theor3r^ and also observed experimentally.'^ 



2. Photo -assisted field emission 



C. Laser-induced thermally-enhanced field 



emission 



Even if Nfiuj < ^, it is possible for electrons to tun- 
nel through the potential barrier if a DC bias is applied. 
The emitted current is described by the Fowler-Nordheim 
equation with an effective work function ^eff = (p — Nfiuo. 
This process is called photo-assisted field emission (Fig. 
[2|b)), and has been observed for laser-illuminated tung- 
sten tips both with continuous illuminatiorP^ ^nd for il- 
lumination with a femtosecond laser. ^ Since the excited 
electron states in metals have short lifetimes, the emis- 
sion process is expected to be prompt, i.e. should fol- 
low the laser intensity profile to first approximation, as 
should multi-photon emission. 



Illumination of the tip by a train of ultrafast pulses 
heats the tip and causes the apex electron temperature 
to evolve periodically in time. Therefore, thermally en- 
hanced field emission can increase drastically at short 
times following laser excitation and before the electron 
gas equilibrates with the lattice. It is natural to divide 
the thermally enhanced current into two components: 
one is the DC offset (DC TFE), and the other is the 
pulsed component (transient TFE). 



1. Transient TFE 



3. Above-threshold photoemission and optical field emission 

As pulse fiuence increases, above-threshold photoemis- 
sion (TV-photon emission with N greater than the min- 
imum number of photons required to escape the metal) 
can occur (Fig. |2|d)). This process is analogous to 
above-threshold ionization of atoms. ^ For sufficiently 
high laser intensities, the lower photon orders may be 
suppressed due to the quiver energy of the electron in 
the laser field. This strong-field effect has recently been 
observed for ultrafast laser-induced emission from golcP^ 
and tungsten tips.^^ 

For sufficiently large peak intensities, the description of 
the process becomes particularly simple: the laser field, 
Fiaser IS stroug cuough that electrons can tunnel from the 
Fermi level during part of the optical cycle, and one can 
substitute F ^ F -\- Fiaser in the Fowler-Nordheim equa- 
tion to estimate the laser-triggered current. This picture 
is applicable when the Keldysh (or adiabaticity) param- 
eter, 7, is less than one.«^ This parameter is defined as 



7 = 



(7) 



The temporal profile of transient TFE depends on the 
electron-phonon coupling time, Tei-ph, and the dynam- 
ics of heat conduction in the tip. The electron gas is 
out of equilibrium with the lattice at timescales shorter 
than Tei-ph, which is typically on the order of picosec- 
onds. Since the electron contribution to the heat capac- 
ity is smaller than that of the lattice, for pulses shorter 
than Tei-ph the electrons temporarily reach high temper- 
atures. The high electron temperature is established on 
a very short (fs) timescale, and its dynamics can then be 
described by the two-temperature model.l^ The electron 
temperature decays by thermal conduction away from 
the apex and by coupling to the phonons; whichever 
timescale is faster will determine the duration of the TFE 
transient. 



2. DC TFE 

Following the transient, the temperature of the tip 
apex is nearly constant for most of the time between laser 
pulses. A temperature offset can arise from the inability 
of the nanostructure to dissipate the deposited energy 
quickly enough. The temperature rise can be estimated 
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as that due to a continuous laser and is approximately 
given b}0^ 



AT ; 



(8) 



where (3 is the fraction of light incident on the tip that is 
absorbed, Pq is the average laser power, k is the thermal 
conductivity of the tip material, wq is the beam waist, 
and is the opening half-angle of the tip, which is conical 
approaching the apex. From AT, the average DC TFE 
current can be estimated from Eq. [l] 



IV. RESULTS AND DISCUSSION 

A. Two-photon emission vs. laser-induced 
thermally-enhanced field emission 

In HfC, the laser- induced emission process exhibits two 
regimes; in this Section we present experimental evidence 
that two-photon emission dominates the current for low 
optical power and bias field, whereas a DC current from 
thermal-assisted field emission (DC TFE) dominates at 
high optical power and bias field. In the range of our 
experimental parameters, the possibility of emission on a 
picosecond time-scale due to transient heating of the elec- 
tron gas (transient TFE) is inconsistent with power and 
field dependence measurements as well as simple pump- 
probe measurements. 

First, the emission regimes are demonstrated as either 
pulsed (at low bias) or DC (at high bias) by comparing 
^laser with as the bias field is varied, as shown in Fig. 
[Sja). The crossover between these regimes occurs around 
1.8 GV/m; above this bias field, /2f accounts for less than 

half of /laser- 

In the low-bias regime, we attribute the pulsed cur- 
rent to A^-photon emision over the potential barrier, with 
N = 2 the dominant contribution. Evidence for two- 
photon emission includes both the bias field and inten- 
sity dependence of /2f (and, equivalently in this regime, 
/laser) which are plotted in Fig. [sj^a) and the inset. The 
inset shows a quadratic intensity dependence over the full 
range of bias fields applied, indicating a two-photon pro- 
cess. Fig. [sja) plots the field dependence for only two 
specific intensities for clarity. The /2f data are fit, using 
a single normalization coefficient, to the model 



/2-ph ^ P''[{Nnuj-c^^^eF/{A^eo)f^^{kBT)\ (9) 

derived from Eq. |6] in the limit that ksT <C {NTiuo — 
(j)). To determine the work function and the factor k 
relating voltage to bias field, the room-temperature I-V 
data (no laser illumination) is first fit to Eq. [sjas shown in 
Fig. [sj^b). The temperature is determined from /laser as 
described in the next paragraph. Then, the two fitting 
parameters of Eq. |9] are N and uj. Choosing N = 2 
results in huj = 1.46 eV, close to the center frequency of 
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FIG. 3. Ultrafast laser-triggered emission from a HfC field 
emission tip. (a) Closed symbols are /laser and open symbols 
are /2f, for P = 3.9 mW (triangles) and 9 mW (circles) laser 
power. The corresponding peak intensity in the focus, ne- 
glecting field enhancement, is 2.1x10^° and 4.9x10^° W/cm^. 
Light solid lines are theory for TFE and dark solid lines are 
fits of Eq. (|9}, describing two-photon emission, to the data. 
Inset: Current vs. average laser power for tip bias of 200-700 
V, increasing in 100 V steps, (b) Cold field emission (/dc at 
300 K) without the incident laser. Theory for 300 K emission 
is fit to the data to determine the work function and emission 
area. Then the family of TFE curves (solid lines, theory) de- 
scribe emission for a tip with elevated temperature, obtained 
using Eq. [l] To replicate the lock-in signal expected if total 
current is dominated by DC TFE at a given temperature, the 
300 K curve is differenced from the TFE curve, to give the 
light solid curves in part (a) of the figure. 



the (broad) laser spectrum. Larger values of N would 
lead to photon energies far outside the spectrum of the 
laser. 

Next, we attribute the high-bias regime dominated by 
DC current to the steady state TFE discussed in sec- 
tion 



III C[ caused by laser heating and the tip's inabil- 



ity to cool to room temperature between incident op- 
tical pulses. Above 1.8 GV/m, /laser is well-described 
by fitting to tempe rature-dependent field emission the- 
ory (Section |III A| ) as shown in Fig. [3]^a). The fit- 
ting is performed as follows: first, the family of TFE 
curves Itfe{T, F) is calibrated as in Fig. [sj^b) by fit- 
ting room-temperature I-V data to the temperature- 
dependent field emission theory described in Section 



III A The curves have been obtained by numerically in- 
tegrating Eq. [l] At a given temperature, the expected 
lock-in signal due to TFE is then calculated by subtract- 
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ing room-temperature current from elevated-temperature 
current. The predicted lock-in signal for different as- 
sumed tip apex temperatures is shown as the light solid 
lines in Fig. [sj^a). The apex temperatures in the experi- 
ment may be inferred to ± 50 K with this method, and 
agree with the prediction of Eq. [S] for our experimental 
parameters, as will be discussed in Section \LV B| 

One may consider whether transient TFE (due to emis- 
sion from the hot electron distribution on a picosecond 
timescale after the laser pulse, as described in section 



IIIC) contributes significantly to current in the pulsed 
regime. In that case, the emission timescale could be on 
the order of the electron-phonon coupling time (~1 ps), 
significantly longer than the ~10 fs expected for two- 
photon emission. However, transient TFE can be ex- 
cluded as the dominant current source for several reasons. 
First, we observe an additive baseline for autocorrelation 
traces measured with the tip as a detectoi!^, as shown in 
Fig. [4[ This is seen for time delays At ^ 50 fs, to avoid 
fringes from the interference of the two pulse trains at 
the tip. Transient TFE should not exhibit such linearity 
when the total energy deposited is doubled. Secondly, the 
expected field and power dependence of transient TFE is 
qualitatively different from that of two-photon emission. 
Current due to transient TFE would scale roughly like 
the TFE curve (Fig. [sj^b)) for the corresponding peak 
electron temperature following the laser pulse, up to a 
normalization factor due to the finite electron-phonon 
coupling timescale. For these curves, the nonlinearity 
of TFE current as a function of laser power changes with 
bias field. This is in stark contrast to the observed 
behavior of /2f , independent of bias field, and can be im- 
mediately visualized by noting that whereas two-photon 
emission curves vs. field are parallel for various laser 
powers on a semi- logarithmic plot (Fig. [sj^a)), the family 
of TFE curves (Fig. ^h)) are not. 

The effect of laser polarization also supports two- 
photon emission in favor of transient TFE. For laser po- 
larization perpendicular to the tip axis, /2f is suppressed 
by a factor of ~25 independent of the bias (Fig. [sja)) rel- 
ative to parallel polarization. This suppression is consis- 
tent with the idea that photons excite electrons in surface 
states, where momentum conservation favors absorption 
of light with polarization vector normal to the surface^^. 
Light polarized parallel to the tip axis satisfies this con- 
dition near the apex, where both DC and optical field 
enhancement are largest. As will be discussed in Sec- 



tion |IVB] the transient temperature rise is expected to be 
larger for parallel polarization, so transient TFE would 
also be suppressed by perpendicular polarization. How- 
ever, the measured suppression is a constant factor inde- 
pendent of bias, which is again inconsistent with thermal 
emission at two different temperatures. 

We also note that the laser heating leading to TFE in 
the high-bias regime is greatest when the beam is aligned 
slightly up the shank^^ because more energy is absorbed. 
In Fig.jsjb), such an alignment change (~2-3 pm) results 
in a much higher temperature. Because the two-photon 
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FIG. 4. /laser as a function of time delay between two identi- 
cal pulse trains incident on the tip. For At ^ 50 fs (well away 
from interference between the two pulses), /laser is indepen- 
dent of the time delay. At A, one of the two pulse trains is 
blocked, showing that the signal decreases to half of its value 
(within ^5%, due to a small power imbalance or difference in 
alignment between the two pulse trains). At B, both beams 
are blocked. The average power is 4 mW per pulse train, and 
the tip bias field is 0.4 GV/m. In this parameter regime, /2f 
and /laser are equivalent. 
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FIG. 5. (a) Polarization strongly influences the two-photon 
emission. Circles (triangles) are for polarization parallel (per- 
pendicular) to the tip axis. Closed (open) symbols are /laser 
(/2f). The dashed line is the calculated TFE current for the 
tip at 1100 K. The ratio of I^i for the two polarization cases 
is constant, (b) Misalignment can excessively heat the tip. 
Circles (blue) are /laser with the laser aligned to the apex of 
the tip, whereas diamonds (red) are /laser measured with the 
laser aligned further up the shank, for the same laser power. 
Light solid lines are DC TFE theory at the labeled temper- 
atures. Alignment up the shank increases the TFE current 
due to the higher apex temperature, but decreases the two- 
photon emission, observable at low fields. The average laser 
power is 9 mW for all of the curves (peak intensity in the 
focus, 4.9 X 10^° W/cm^) 



7 



multi-photon 



thermal 



FIG. 6. Emission pattern under laser illumination for two 
different DC bias fields applied to the tip. At lower bias fields 
(left), multi-photon emission dominates while at higher bias 
fields (right) thermally-enhanced field emission takes over. 
Cold field emission with no illumination gives a pattern iden- 
tical to the TFE whereas the multiphoton emission is more 
diffuse and emits from slightly different crystal planes. 




FIG. 7. Images of the HfC tip taken with a scanning electron 
microscope. The scale bars, from left to right, are 0.5, 2, 
and 10 pm. The viewing perspective is such that the tip is 
directed out of the plane by ^ 45deg. 



emission (the dominant contribution to /laser for low bias 
fields) depends mostly on the intensity of light at the 
apex, it decreases when the beam is aligned up the shank 
and the two curves in Fig. [5|b) cross each other near 0.5 
GV/m. All other data presented in this work is measured 
with the beam aligned to the apex by maximizing /2f- 

Besides the identification of the two emission processes 
via I-V measurements, we note that the emission pattern 
is significantly different for the two processes (Fig. [6]). 
When TFE dominates, the pattern is essentially iden- 
tical to the DC emission pattern without laser illumina- 
tion. However, when two-photon emission dominates, the 
pattern is more diffuse and comes from different crystal 
planes. We have observed this for each of several differ- 
ent 310-oriented HfC tips, allowing for the rapid identi- 
fication of the dominant emission mechanism. We note 
that this differs from the case of photofield emission from 
Ill-oriented tungsteiP, where the prompt laser-induced 
emission has the same emission pattern as DC emission. 



B. Model of laser heating and temperature 
dynamics 

A finite element analysis model^^ of the evolution of 
the temperature at the tip apex following laser illumina- 
tion is developed. The model allows us to predict how 
transient TFE scales with laser power and bias field in 
comparison wit h the prompt two-photon emission of in- 
terest (see Sec. JVC). The first stage of the model is a 
solution for the steady-state temperature distribution in 
the tip under conditions of equilibrium between the elec- 
tron gas and the lattice. This is used as the initial state 



TABLE I. Materials properties used for modeling. 



[37W391 



T (K) HfC 



W 



Au 



C (10^ J/(m'' K)) 
Cei (10^ J/(m3 K)) 

[42lf43l 

K (W/(m K)) 

A.ei,o (W/(m K)f^ 
^ei-ph (10^^ W/(m3 K)pl 



300 
1000 



2.29 
3.48 



2.55 
2.89 



2.48 
2.88 



P (II polarization) 
/3 (± polarization) 



300 3.8r 
1000 8.4^" 


3.46 2.02 
9.88 6.73 


300 20.7 
1000 27.0 


171 317 
117 270 


b 


ISr 318 


6.85| 


^ 19.1 2.61 


Re[] 9.4^ 
Im[] -8.48 


T 5.2 -24 
^ -19.4 -1.5 


0.8 
0.95 


0.6 0.037 
0.7 0.047 



^ HfC Cei(T) determined from the density of states in Ref. 1401 
^ For HfC, ^Cei and ^Cph were taken to be 12 W/(m K), half of k 

at 630 K, the approximate initial temperature (see text). 
^ W ^Cel cit 400 K determined by fitting electrical resistivity and 

thermal conductivity data in the 300-900 K range (Ref. 43 ) 

using the Wiedemann-Franz law. 
d Example value used in HfC TTM model. 

^ Determined by Drude model fit to reflectivity curves in Ref. 1451 
using two Lorentz oscillator terms. 



of the tip before the arrival of the next laser pulse in the 
second stage of the model. The second stage solves for 
the evolution of Tqi in the first ns after the arrival of the 
laser pulse. Representative material parameters used in 
the models are given in Table |lj 

The tip geometry is based on SEM characterization of 
the tip used in the experiments; for the 120-nm radius tip 
imaged in Fig. [7| the upper 1.8 microns of the model are 
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FIG. 8. Modeling the laser- induced heating. An incident 1 
V/m monochromatic plane wave is assumed, and the har- 
monic response of the electric fields is solved with finite ele- 
ment analysis software, (a) Ex at the particular phase when 
it reaches its maximum for an x-polarized incident wave, (b) 
Cycle-averaged magnitude of the electric field for incident x- 
polarized light, (c), (d) Average power density heating the 
tip for X- and y-polarized incident light, respectively. Electric 
field units are V/m and power density units are 10^ W/m^. 
All plots are for the cross section given by the symmetry plane 
(y = 0) of the model. 



shown in Fig. [S] The 2.5 degree half- angle cone extends 
to 16 pm below the apex before slowly flaring out to a 
500 pm diameter cylinder held at 300 K. A hybrid model 
was used for the spatial distribution of heating in the tip 
in order to make feasible the computation of heating near 
the apex, where optical field enhancement^^ has com- 
plicated effects. Maxwell's equations are solved near the 
tip apex, while a simple calculation of the energy inter- 
sected by the tip's silhouette is used as the heat source 
further up the shank, as described below. 

Fig. |8] shows how the detailed heating distribution 
is modeled near the tip apex for an incident 800-nm 
plane wave. Fig. [sj^a) depicts the axial component of 
the electric field at the phase when it reaches its max- 
imum strength (field enhancement 1.9), and Fig. [sj^b) 
shows the cycle- averaged amplitude of the electric field, 
both for x-polarized incident light. The plane wave of 
amplitude 1 V/m is incident from the left side of the fig- 
ure, and the highest time-averaged fields occur toward 
the shadow side of the apex in agreement with Ref. 48, 
The accuracy of the model was checked by replicating 
the models in Ref. 48 and by benchmarking against Mie 
scattering calculations for a 100-nm gold sphere. 

The field solutions are used to calculate the power 
density of resistive heating inside the tip, Pmodei(^, 
as shown in Fig. [8jc,d) for x- and ^/-polarized incident 
light, respectively. We approximate the resistive heating 
caused by the broad spectrum of the laser by the heat- 
ing due to 800 nm light. The integrated resistive heating 
changes by less than 10% if the problem is solved for ei- 
ther 680 nm or 940 nm incident light, and the distribution 
does not noticeably change. 

Although the electromagnetic model is solved using an 
incident plane wave, the experiment uses a focused beam. 



We approximate the heat source in the apex region as 

Papex(^, Z) = Pmodelix, y, z) g-^C^/^o)' (^Q) 

since the beam waist is much larger than the size scale 
of the tip geometry. Further up the shank, the time- 
averaged volumetric heat source Pshank(^) is taken to be 



Pshank(^) 



:wQr{x) 



-2{x/woY 



(11) 



where r{x) is the radius of the tip shank at distance 
X from the apex and Pq is the laser's average power. 
The constant factor /3 takes into account the interrelated 
effects of field enhancement and reflectivity of the tip 
material, and is determined by requiring that Pshank(^) 
smoothly match the numerical solution (integrated over 

z) near the apex. Although the expression for Pshank 
ignores the finite penetration depth of heating due to 
skin effects, details of the shank heating distribution af- 
fect neither the equilibrium temperature distribution nor 
the short-timescale dynamics of the temperature near the 
apex, which are the quantities of interest. 

To model the steady-state temperature rise of the tip 
during laser heating, we use a continuous (i.e. not pulsed) 
heat source, 

p(a:,2/,.) = ^-^'^("'^'")' -<l-5l^"^ (12) 

[Pshank(^), X > 1.5 pm 

for t > 0. Temperature-dependent values of the thermal 
conductivity and heat capacity of HfC were taken from 
the literaturel ^^ l ^^ l The model predicts a DC temperature 
rise of 262 K and 567 K for average laser power of 3.9 
mW and 9 mW, respectively. These values are in good 
agreement with the observed values of 225±50 K and 
550±50 K (see Fig. [3|. The simple model of Eq. [s] 
predicts apex temperature increases of 300 K and 700 K 
for 3.9 and 9 mW, with = 2.5 degrees and (3 = 0.95 
(note that Eq. [s] assumes a simpler tip geometry). The 
time-scale for reaching the equilibrium temperature after 
the heat source is turned on is tens of microseconds. 

The second stage of the model investigates the effect 
of a single laser pulse once the steady-state temperature 
profile has been established. The electron and lattice 
temperatures (Tei and Tph) are described by the coupled 
equations 

= V • (A>:elVTei) - ^el-ph(Tel - Tph) +p(:^,^, z)f{t) 

(13) 

= V • (A>:phVTph) + ^ei-ph(Tei - Tph), (14) 



a 



where ^ei-ph is the electron-phonon coupling constant. Tq\ 
and Tph are functions of position and time, and Cg, ph and 
A^e, ph may be functions of position and time via their de- 
pendence on Tel and Tph- The temporal profile of the 
heat source, /(t), is on at a constant value for < t < 10 
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FIG. 9. Results of the TTM simulation for HfC (details of 
the model are provided in the text). Upper plot: Temper- 
ature at the tip apex as a function of time. Black (darker) 
and red (lighter) curves are for x- and ^/-polarization, respec- 
tively. Thick solid curves are Tei and thin dashed curves are 
Tph. Thin dotted curves show Tei without electron-phonon 
coupling for comparison. The inset is rescaled to show Tph 
more clearly. Lower plots: the Tgi distribution (?/ = cross- 
section) at several times after a laser pulse. For clarity Tscaie, 
the maximum temperature of the color scale, is 1250 K for 
the 10 fs and 1 ps frames, and 640 K otherwise. 



fs and zero otherwise. It is normalized to represent the 
heating from a single pulse in the 150 MHz pulse train for 
4.6 mW average power. The initial condition for Tei and 
Tph is the temperature profile obtained for steady-state 
heating with 4.6 mW average power. The electron ther- 
mal conductivity for HfC is unreported in the literature 
but is expected to be similar to that of TiC, approxi- 
mately half the total thermal conductivity.^^ Since the 
electron-phonon coupling time is unknown, two solutions 
are shown in Fig. [9j the electron temperature with no 
electron-phonon coupling (^ei-ph = 0; thin dotted curves) 
and solutions for both Tei and Tph (thick curves and thin 
dashed curves, respectively) assuming ^ei-ph = 6.85 x 10^^ 
W/(K m^), independent of temperature, which corre- 
sponds to Tei-ph ~ C'ei/^ei-ph ~ 1 ps. The solutiou with 
no electron-phonon coupling is an upper limit to the elec- 
tron temperature. The apex temperature in this solution 
decays via conduction on a timescale of hundreds of ps; 
normal values of electron-phonon coupling times are on 
the order of picoseconds, which implies that (barring an 
anomalously long Tei-ph) the duration of TFE from the 
HfC tip following a laser pulse is set principally by the 
electron-phonon coupling time. The electron tempera- 
ture distribution near the tip apex for the HfC model 
with electron-phonon coupling is also shown in Fig. |9] 



HL 



FIG. 10. (Top row) Modeled, cycle- averaged electric field 
strength near the tip apex for a 1 V/m, 800 nm, x-polarized 
incident plane wave, ^scaie is 2.1, 2.2, and 2.4 V/m for HfC, 
W, and Au, respectively. (Bottom row) Energy density de- 
posited by a single laser pulse in the 4.6 mW, 150 MHz pulse 
train when focused to a 2.7 pm waist. The color scale ranges 
from to tscaie = 5.0, 8.5, and 0.8 x 10^ J/m^ for the three 
materials, respectively. 



Note the change in the temperature scale between the 
1 ps and 3 ps plots as the electron temperature decreases 
due to electron-phonon coupling without significant spa- 
tial redistribution of thermal energy from conduction. 
The model can be used to compare HfC, tungsten, and 



gold as tip materials. Fig. [TO] shows the time-averaged 
magnitude of the electric field near the apex and the heat 
deposited by a single pulse in the 4.6 mW pulse train, for 
x-polarized incident light. Although the field enhance- 
ment outside the tip is only affected at the 10% level by 
the various materials' permittivities, the heat density de- 
posited at the apex varies by an order of magnitude and 
the total heat deposited in the gold tip is 26 times less 
than for HfC. However, the smaller skin depth, but sim- 
ilar reflectivity, of tungsten as compared to HfC leads to 
a 70% higher heat density impulsively deposited in the 
electron gas at the apex. 



Fig. [TT] shows the modeled apex electron temperature 
due to transient heating from a laser pulse for tungsten, 
and gold with the same experimental parameters as Fig. 
|9]and material parameters in Table [T[ The electron ther- 
mal conductivity for W and Au as a function of electron 
and lattice temperature is expanded aJ^ 



/^el(^eb^ph) = ^el,0 



^ph 



(15) 



Conduction of the lattice is insignificant and therefore 
neglected (^^ph = 0). The values of ^ei-ph (see Table |l| 
are taken from Ref. ^2 averaged over the temperature 
range of the simulation. Several qualitative trends are 
captured in Fig. 
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First, due to gold's high reflectiv- 
ity and both gold and tungsten's high thermal conduc- 
tivities, the initial steady-state elevated apex tempera- 
ture is very low compared to that of HfC. Therefore DC 
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Time [s] Pulse Energy (nJ) 



FIG. 11. Results of the TTM simulation for W and Au. The 
curves are plotted in the same way as in Fig. [5] 

TFE is essentially negligible in tungsten and gold for the 
range of average optical power used in high repetition 
rate experiments.^^ Second, x-polarization heats the 
tip apex more than ^/-polarization for all three materials, 
as observed by the initial temperature rise for the single 
pulse. For ^/-polarization in HfC, the heat is primarily 
deposited slightly away from the apex (see Fig. jsjd)), so 
that, without electron-phonon coupling, the apex tem- 
perature continues to rise after the optical pulse due to 
thermal conduction. Third, for gold and tungsten, the 
apex cools via diffusion much more rapidly than HfC in 
the absence of electron-phonon coupling. This is due to 
both the higher electronic thermal conductivities of these 
two materials and the fact that the heating is initially 
more narrowly confined (smaller optical skin depth). In 
gold, because the mean free path of the excited electrons 
is larger than the optical skin depth, ballistic tr anspo rt 
would reduce the apex temperature even further.^^^Ml 



C. Bounds on transient thermally-enhanced field 
emission 

The fractional contribution of transient TFE to the 
pulsed current (measured by l2^)^ as well as its scaling 
with laser power, should be compared with that of multi- 
photon emission in order to understand the regimes in 
which the tip can be operated as a fs-timescale emit- 
ter. We note that the field dependence of /2f (Figs. Isl 
[5| is accurately fit by two-photon emission using Eq. ^ 
throughout the parameter regime accessible in our exper- 
iment (limited to P ^ 11 mW, above which the degree of 
DC TFE induced by small mis-alignments of the beam is 
much larger than the pulsed current). This suggests that 



FIG. 12. Upper bound to peak temperature at the tip apex 
immediately after heating by a single pulse in a 150 MHz 
pulse train of 10 mW average power focused to 2.7 pm, as 
a function of assumed Te-ph- The inset provides an exam- 
ple calculation: Open circles are /2f and closed circles are 
/laser. The two soUd curves represent transient TFE at dif- 
ferent temperatures, scaled by Te-ph so as not to exceed /2f. 
The steeper curve is 1200 K emission for 257 ps, and the more 
shallow curve is 2000 K emission for only 34 fs. (b) Scaling of 
two-photon emission (solid lines) and the estimated transient 
TFE (dashed curves) for a single pulse incident on a room- 
temperature tip. The lower (blue) and upper (red) curve 
in each pair are for a 0.5 and 1.5 GV/m bias field, respec- 
tively, the heat deposited per nJ pulse energy was taken from 
the FEM model, and 1 ps electron-phonon coupling time has 
been assumed. The largest pulse energy shown corresponds 
to 7 1.5, ignoring the Schottky barrier reduction. 

transient TFE contributes negligibly to the total pulsed 
current. 

We can place the most stringent bound on the 
contribution of transient TFE using experimental 
data for 10 mW (the highest power measured) with 
perpendicularly-polarized light (suppressing two-photon 
emission by a factor of ~ 30). This requires an es- 
timate of the field dependence of the transient TFE, 
as follows: First, the measured DC thermal emission 
at room temperature is used to calibrate the family of 
TFE curves Itfe{T^ F) for this particular data set, as 
in Fig. |3|b). Next, a conservative upper limit to the 
transient TFE as a function of bias field is estimated as 
a square current pulse lasting for a time Tei-ph, due to 
the peak electron temperature before coupling with the 
phonons. Therefore, the upper limit to transient TFE is 

-^transient 

(F) < /TFE(^max,^)Tel-ph/rep < hi{F) , where 

average current has been used for comparison with exper- 
imental I-V data. Given an assumed value of Tei.ph, ^max 
is the maximum electron temperature at the apex such 
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that /transient Calculated in this way does not exceed the 
measured /2f at any bias field. Fig. [l2|a) shows T^ax as 
a function of Tei-ph- For rei_ph = 1 ps, T^ax ~ 1690 K and, 
for comparison, the peak apex temperature predicted by 
the modeling described in Section \LVB is well below this 



at 1490 K.l^The inset to Fig. 
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^a) shows two examples 

of /transient (/") choSCU tO illustratC hoW T^ax dcpcuds OU 

the value assumed for Tei-ph* a relatively high electron 
temperature is only consistent with experimental data 
(/2f) assuming a short Tei-ph, whereas a relatively cool 
temperature could emit for a much longer coupling time. 
Finally, both the two-photon and transient TFE contri- 
butions to the pulsed current must be scaled to account 
for using parallel polarization in any experiment designed 
to maximize the flux from fs electron pulses. The ~30- 
fold increase in two-photon emission was directly mea- 
sured, as in Fig. [5|b), whereas modeling predicts that 
parallel-polarized light deposits ~65% more energy den- 
sity in the apex than perpendicularly-polarized light. 
Taking both scalings into account, over the useful range 
of bias field {F > 0.5 GV/m where the two-photon emis- 
sion turns on strongly, and F < 1.4 GV/m, where DC 
TFE becomes important), we estimate that transient 
TFE could have contributed no more than 5% to 25% 
of the two-photon signal. We emphasize that this is a 
conservative upper bound, since it corresponds to the 
situation shown in the inset of Fig. [l2|a), where the 
transient TFE would account for all of /2f precisely at 
the first or last data point measured. 

Similar ideas can be used to estimate how the transient 
TFE scales with laser pulse fluence in comparison with 
the quadratic two-photon emission scaling (Fig. p!2|b)). 
In this experiment, the high repetition rate results in ap- 
preciable DC heating so that DC TFE dominates the cur- 
rent before there is appreciable transient TFE. However, 
in an experiment requiring large peak current densities, 
the laser pulse energy could be increased while decreas- 
ing the repetition rate to avoid DC heating. In the limit 
of low repetition rates, the steady-state temperature of 
the tip is assumed to be 300 K. The initial temperature 
rise of the electron gas is estimated using an approximate 
energy density, Q = PqT^ with T the laser pulse fluence, 
and /3q = 220 J/m^ per nJ/cm^ determined from the 
model of section HV B[ The maximum transient electron 
temperature at the apex of the tip, Ttrans is determined 
by solving 



Q 



r 



C,i{T)dT. 



(16) 



As before, the normalization of the two-photon and 
TFE curves is determined from experimental data for 
our tip. An upper limit to the average current due to 
transient TFE is calculated as a square pulse of dura- 
tion Tei-ph, emitted from the apex with electron tem- 
perature Ttrans- Explicitly, Fig. [l2jb) shows /trans = 
/tFE (/trans /^j^el-ph/rep aS a fuUCtioU of pulsC fluCUCe, 

for two values of bias field. We find that transient TFE 
is negligible over a wide range of pulse fluences, though 



it may become important as the power increases to the 
intermediate regime (7 ~ 1).-^ 



V. CONCLUSION 

We have measured two-photon emission from a 
hafnium carbide tip induced by sub- 10 fs laser pulses. 
HfC's low work function and demonstrated stability for 
high current density emission make it an excellent can- 
didate tip material for ultrafast laser-triggered electron 
emission applications, particularly those with constraints 
on available peak laser intensity. At a repetition rate of 
150 MHz, the maximum current achieved is 5 electrons 
per pulse, limited by the onset of thermally-enhanced 
field emission due to an increase in the average tempera- 
ture of the tip. We estimate that for light polarized par- 
allel to the axis of the tip under the conditions studied, 
transient TFE contributes negligibly to the total current. 
We present a simple framework for understanding the 
relative scaling of multiphoton emission and thermally- 
enhanced field emission with changing laser parameters 
and bias field. 

There are two complementary regimes for time- 
resolved pump probe measurements involving electrons. 
One end of the spectrum uses large numbers of elec- 
trons in a pulse in order to achieve single-shot measure- 
ments. At the opposite extreme, space-charge broad- 
ening may be completely eliminated by operating in 
a single-electron-per-pulse regime, with high repetition 
rates to reduce acquisition time. Both of these regimes 
may be affected by the thermal effects studied here, at 
least in principle: at high repetition rates, the DC tem- 
perature rise of the tip can result in significant DC TFE, 
while at very high pulse fluence, one should consider 
the possibility of large temperature rises for t < r^-ph 
and accompanying transient TFE. The transient effect 
is more insidious since the resulting picosecond-timescale 
pulse would still be below the time resolution of stan- 
dard electronic detection. However, we show that it is 
not a significant limiting factor under the conditions of 
our high repetition rate experiment, even in HfC with its 
high degree of laser heating and relatively poor thermal 
conduction. Comparison of tip materials shows that the 
favorable thermal properties of gold and tungsten greatly 
reduce both the steady-state temperature rise of the tip 
and the timescale of transient emission. We have shown 
that HfC can be operated with relatively little optical 
power as a high repetition rate ultrafast point source of 
electrons in the technologically attractive regime of single 
electrons per pulse. 
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